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Abstract 

We review recent theoretical developments in understanding the many-body per- 
turbative QCD dynamics of strongly-interacting hard probes in dense nuclear mat- 
ter. The relation between initial- and final-state parton scattering in nuclei and 
the quark-gluon plasma and experimental measurements of the Cronin effect, nu- 
clear shadowing, jet quenching, modified di-hadron correlations and forward rapidity 
particle suppression in p+A and A+A reactions is clarified. Our approach empha- 
sizes the process dependence of nuclear effects and outlines the techniques for their 
dynamical calculation and incorporation in simulations that build upon collinear 
factorization in QCD. 
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1 Introduction 



In high energy p+p collisions jet and particle production is well described in 
the framework of the collinear PQCD factorization approach. However, when 
a partonic process is embedded in nuclear matter, the cross sections are modi- 
fied by the soft scatterings that precede and/or follow the hard interaction [1]. 
One can gain insight in the various aspects of perturbative many-body QCD 
dynamics by following the history of a parton from early times, t — > — oo, in the 
parton distribution function to its asymptotic state, t — * oo, in the wave func- 
tion of the observed hadron. In our approach, the strength of jet interactions 
is controlled by microscopic and non-perturbative parameters characterizing 
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Type of scattering 


Pt dependence of the nuclear effect 


Elastic (incoherent) 
Inelastic (final-state) 

Inelastic (initial-state) 
Coherent (t—, u— channel) 
Coherent (s— channel) 


IIP: Cronin effect; I2P: di-jet acoplanarity 

IIP: quenching at all pt', I2P: quenching/enhancement 

at high/low pt 

IIP: suppression at all pr] I2P: same as IIP 
IIP: suppression at low p?; I2P: same as IIP 
IIP: enhancement at low pr; I2P: same as IIP 



Table 1 

Effect of elastic, inelastic and coherent multiple scattering on the transverse mo- 
mentum dependence of single (IIP) and double inclusive (I2P) hadron production 
in the perturbative regime. 

the medium, such as temperature T, energy density e, the scale of high twist 
corrections £ 2 , typical momentum transfer squared /z 2 , and mean free path A. 
Results are summarized [1] in Table 1. 



2 Cronin effect and dynamical shadowing 



We first consider partonic 2 — > 2 processes, which we call elastic as opposed to 
inelastic 2 — > 2 + n processes associated with gluon bremsstrahlung. Multiple 
interactions of the incoming quarks and gluons lead to transverse momentum 
diffusion [2]. For initial parton flux d 2 N l (k) = 5 2 (k)d 2 k the distribution of 
partons at the hard collision vertex reads: 
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Such broadening can be implemented in the PQCD formalism as follows: 

(^) = (kLc) + 2x l /i 2 e, < = 1,2, (2) 



and leads to y^-dependent enhancement of the hadron cross sections at inter- 
mediate pt ~ few GeV. The left panel of Fig. 1 shows theoretical calculations 
at y = compared to p+W/p+Be data and predictions for d+Au collisions 
at RHIC. We used 2/j, 2 /\ q w 0.12 GeV 2 /fm, 2/j, 2 /\ g w 0.28 GeV 2 /fm. For a 
survey of models of the Cronin effect see [3]. 
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Pt [GeV] 

Fig. 1. Left panel: the ratio of ^4-scaled p+W/p+Be data on 7r + ,7r~ production 
at = 27.4,38.8 GeV. Right panel: all-twist resummed F2(A) / F2(D) calculation. 
The band corresponds to the choice ^ 2 = 0.09 — 0.12 GeV 2 . 

A distinctly different effect arises from coherent final-state interactions in large 
nuclei [4]. When the light cone momentum fraction xb < 0.1 the longitudinal 
extent of the exchange virtual meson l c = 1/2x^771^ exceeds the nucleon size 
tq and the struck parton scatters coherently on several nucleons. The high 
twist scale per nucleon is defined as: 



e = 5^p( P | F 2 (\j) W = 3 -^fl l xG(x) . (3) 

Nuclear shadowing can be cleanly studied only in deep inelastic scattering 
(DIS). Resummation of all nuclear-enhanced power corrections to the inclusive 
cross sections generates dynamical mass for the final-state parton m\ yn = 
£2^1/3 Q n exam pi e f t ne transverse structure function we have: 



Cft* A Fr (* + O 2 ) = F™ (« + . (4) 

The left panel of Fig. 1 shows good description of the world's data on DIS 
on nuclei [4] with £ 2 = 0.09 — 0.12 GeV 2 and Q 2 > m 2 N . Similar results 
were derived for p+A reactions [5] and our calculations have been generalized 
to include the case of heavy quarks [1]. For a survey of models of nuclear 
shadowing see [3]. 
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Initial-state multiple scattering in the coherent regime, in contrast to final- 
state, leads to enhancement of the cross sections and Cronin effect [1]. Such 
distinctly different behavior points to the process dependence of nuclear effects 
and the need for their dynamical evaluation. 



3 Initial- and final-state energy loss 



The acceleration of the incoming and outgoing partons by the soft interactions 
in the dense nuclear matter results in medium-induced gluon bremsstrahlung [7] 
Given its quantitative importance for the description of the large jet attenu- 
ation in a quark-gluon plasma (QGP) of initial temperature T ~ 400 MeV 
and energy density eo ~ 20 GeV/fm 3 , recent perturbative calculations have 
been focused on final-state radiative energy loss. The goal of such studies is 
to describe the large five-fold suppression of high-px hadron production in 
central Au+Au reactions [8,9]. 

We here rely on the Gyulassy-Levai-Vitev (GLV) approach for describing the 
multiple interactions of the propagating jet +gluon system [10,11] since it was 
specifically developed to describe the QCD dynamics of heavy ion (p+A and 
A+A) reactions. It is well equipped to account for some of the most important 
characteristics of nuclear collisions: L/X g < few <C oo and Ej Ct <C oo. To first 
order in the correlation between multiple scattering centers the differential 
bremsstrahlung spectrum is given by [10]: 
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(5) 



Here q drr . is possible ^ directed momentum transfer, for example transverse 
flow. Clearly, this effect amounts to a change of the reference frame for the 
radiative spectrum but neither distorts it nor changes the amount of lost 
energy AE. 

Integrating over the gluon bremsstrahlung intensity in the limit 2E/ ji 2 L 1 
we find: 



. _ C R a s /2 2 T2l 2E . 9tt«3 1 dN9 2E T 

for static and (1+1)D Bj or ken-expanding plasmas, respectively [10,11]. Tests 
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Fig. 2. Left panel: comparison of the analytic Raa(pt) to PHENIX and STAR data 
on the A^part-dependent quenching. Right panel: nuclear modification factor Raa for 
Au+Au collisions at yfs = 200 GeV versus px and centrality. Similar calculation 
for Cu+Cu at y/a = 200 GeV at RHIC. 

of parton energy loss at RHIC and the LHC can be carried out via the cen- 
trality dependence of jet quenching. With dN 9 /dy oc iVp ar t, L oc and 
A± oc iVp art , the GLV approach gives a definite prediction for the system size 
dependence of energy loss and jet quenching [12]: 

e = AE/E oc NlH , In R AA « -kN^ . (7) 



Our analytic result for the system size dependence of Raa{pt) is shown in the 
left panel of Fig. 2. It is fixed by the magnitude of the suppression established 
in central Au+Au collisions. From this analysis we expect a factor ~ 2 sup- 
pression in central Cu+Cu collisions. Comparison to the PHENIX > 7 GeV 
data [8] and STAR p? > 6 GeV data [9] is also shown. 

To verify the universal suppression pattern, Eq. (7), we carry out numerical 
simulations of Raa versus pt and centrality with QGP properties constrained 
by the entropy density, or equivalently dN 9 /dy and tq at RHIC. The modi- 
fication of inclusive hadron production from final state radiative energy loss 
can be represented as [12]: 

1 — 2 1 

D h/e (z)=> J deP(e) T ^-D h/e (j^+ J de ^(e) -D h/g g ) . (8) 

o z 
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Fig. 3. Left panel: large angle distribution of the radiative gluons relative to th jet 
axis. Right panel: jet softening - redistribution of the energy from high px to low 
Pt particles [14]. Top and bottom panels illustrate central and peripheral collisions, 
respectively. 

In Eq. (8) the first term is associated with the quenching of leading hadrons 
and -P(e) is the probabilistic distribution for the fractional energy loss e = 
AE/E = Y^=i e ii e i = Ui/E, due to multiple gluon emission. The right panel of 
Fig. 2 shows the jet quenching results versus pt for three different centralities 
in Au+Au and Cu+Cu collisions at RHIC [12]. Data is from PHENIX [8]. 

With fractional energy loss as large as ~ 25% for quarks and ~ 50% for 
gluons [12], the quenching of hard jets should be correlated with measurable 
enhanced production of low-pr hadrons [13]. The second term in Eq. (8) re- 
flects the gluon feedback contribution, which for single inclusive measurements 
becomes important only at the LHC. For away-side di-hadron correlations this 
effect is already large at RHIC [14]. If follows from Eq. (5) that the angular 
behavior of the medium-induced gluon bremsstrahlung is significantly broader 
when compared to vacuum radiation. Its most important feature is the cancel- 
lation for the collinear k — > qai r . radiation. The left panel of Fig. 3 shows the 
angular distribution of uj = 3 GeV radiative gluon off a quark jet for L = 5 fm, 
li — l GeV, \ g = 1 fm [14]. The remnant of the leading parton is not shown. 

Depletion of hadrons from the quenching of the parent parton alone leads 
to a large suppression of the double inclusive cross section with weak p T , 2 
dependence, shown in the right panel of Fig. 3. Hadronic feedback from the 
medium-induced gluon radiation, however, completely changes the nuclear 
di-hadron modification factor R AA . It now shows a clear transition from a 
quenching of the away-side jet at high transverse momenta to enhancement 
at pt 2 < 2 GeV [14] compatible with STAR measurements [13]. 

The physics process that has so far not been addressed is initial state energy 
loss in cold nuclear matter [1,15]. To first order in opacity: 
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Fig. 4. Left panel: calculated suppression of 7r° production at y = 4 for d+Au 
collisions at ^/snn = 200 GeV at RHIC. Right panel: nuclear modification of single 
inclusive 7r° (left) and D° + D + mesons (right) in d+Au collisions at RHIC versus 
rapidity and centrality. 



dN„ 



10 



Crol s f dAz 



du ii 2 k 7r 2 



KM 



J 7r(q 



X 



q- 



2 + ^{z)f 
2 kq 



X 



q 



k 2 (k 



k 2 (k- q) 2 



COS 



'k 2 A/ 
2uj 



(9) 



From Eq. (9) with the requirement udN g / 'du d 2 k > we can evaluate the 
energy loss and implement it in the PQCD hadron production formalism: 



/aE L ( x \ 

e = -gT OC - = KA 1/3 , <j) q , g/N {x, Q 2 ) -> 4> q ,g/N [Y^- e > Q 2 ) ■ ( 10 ) 

As shown in the left panel of Fig. 4, shadowing calculations can only account 
for ~ 50% of the hadron suppression at forward rapidities p+A reactions 
at RHIC [1,5]. Taking into account initial-state energy loss, we find good 
description of the y = 4 Raa(pt) for n°s by STAR [16]. The right panel of 
Fig. 4 shows the same perturbative calculation for light hadrons and heavy 
mesons [1] at intermediate rapidities y ~ 2 and ^ = 200 GeV d+Au collisions 
at RHIC compared to PHENIX data [17]. 



4 Conclusions 



Dynamical calculations of initial- and final-state elastic [2], inelastic [10,11] 
and coherent [4,5] multiple parton scattering have been incorporated in the 
framework of the perturbative QCD factorization approach to describe high 
Pt and Et observables in collisions of heavy nuclei. They naturally account 
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for the process dependence of nuclear effects [1] and provide adequate uni- 
fied description of seemingly disparate phenomena, such as the Cronin ef- 
fect [3] and nuclear shadowing [6] in p+A / t + A reactions and high-p T jet 
quenching [7,12] and low pt large-angle di-hadron enhancement [14] in A+A 
reactions. In cold nuclear matter, parton diffusion, nuclear enhanced power 
corrections and initial state energy loss are compatible with similar momen- 
tum transfer scales: ^A 1 ^ 3 « 2/i 2 L/A. In the QGP, jet quenching results are 

2/3 

in agreement with the predicted energy and centrality oc iV^t dependence of 
the non-Abelian parton energy loss. 
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